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SUMMARY 


Experiments have been conducted to determine average skin-friction 
coefficients in the absence of heat transfer for completely turbulent 
flow along the cylindrical portion of cone -cylinder bodies of revolution 
having over -all fineness ratios of 10, 15, and 25. The friction data 
were obtained by directly measuring forces. Numerous boundary-layer sur- 
veys were made to enable all data to be based on an effective starting 
position of the turbulent flow. Mach numbers of 0.5^ 0.8, 2.0, 2.5^ 2.9^ 
3 .^, and 3 *6, and Reynolds numbers between h million and 32 million were 
investigated. At a Mach number of 2.0, data were obtained for different 
pressure distributions by distorting the flexible -plate walls of the wind 
tunnel . 

The results show no significant effect on average skin-friction coef- 
ficient of small changes in pressure gradient. At both subsonic and 
supersonic velocities, the skin-friction coefficient depends only to a 
small extent on cylinder length -diameter ratio. For each length -diameter 
ratio, however, the effect of Mach number is large, amounting to approxi- 
mately a 50 -percent reduction in skin-friction coefficient as the Mach 
number is increased from 0 to 4. This effect of Mach number does not 
depend significantly on Reynolds number or cylinder fineness ratio. 

Comparison with similax experiments on flat plates, primarily those 
of Coles, shows good agreement as to the effect of Mach number on skin 
friction, although the velocity profiles measured on cylinders are 
appreciably different from those on a flat plate. Determinations of 
friction by the momentum method differed considerably from the more reli- 
able direct-force measurements. 


INTRODUCTION 


One of the Interesting aspects of average skin-friction measurements 
for turbulent boundary layers is the unusually extensive history of exper- 
iments conducted in this field. A comprehensive survey of the numerous 
experiments may be found, for example, in the often-cited paper of 
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Schoenherr (ref. l) which covered measurements prior to 1932. More recent 
surveys by Locke (ref. 2) and Hughes (ref. 3) have covered measurements 
made in the subsequent 20 years. An accurate knowledge of skin friction 
for turbulent incompressible flow has, of course, long been of concern 
in the field of naval architecture, with the result that measurements of 
surface skin friction date at least as far back as some experiments of 
Beaufoy conducted near the end of the l8th century. When these early 
data are reduced to dimensionless coefficient form and plotted as a 
function of Reynolds number, the resulting agreement with recent measure- 
ments is surprisingly good. (See refs. 1 and 2.) 

Theoretical investigations of skin friction for turbulent flow, on 
the other hand, are of comparatively recent origin. It is well known that 
the semiempirical mixing-length analyses of von Karmm and Prandtl, 
advanced approximately 20 years ago, yield fairly accurate predictions of 
turbulent friction for a flat plate in low-speed flow. Interestingly 
enough, though, subsequent detailed hot-wire measurements, such as those 
of Liepmann and Laufer (ref. h) and Townsend (ref. 5)^ have indicated the 
fundamental mixing-length hypothesis to be incorrect. The good agreement 
between predicted skin friction and experiment can be attributed in con- 
siderable part to the fact that mixing-length analyses yield arbitrary 
constants of integration which are free to be adjusted by comparison with 
experiment. With constants empirically determined from friction measure- 
ments, the corresponding velocity profiles predicted by mixing-length 
theory using these same constants do not agree well with measurements 
except over a limited portion of the boundary layer. This was noted by 
Gurjienko (ref. 6) and is apparent in almost all measurements dating from 
those of Schultz -Grunow (ref. 7) up to the most recent and detailed 
measurements of Klebanoff (ref. 8) . 

The most satisfactory analytical treatment of the low-speed turbulent 
boundary layer on a flat plate does not utilize mixing-length suppositions, 
but combines instead two empirical laws concerning the similarity of 
velocity profiles: the so-called "wall*' law which was formulated by 

Prandtl, as noted in reference 9 , and the so-called ’Velocity-decrement" 
law which was first observed for pipe flow by Darcy almost 100 years ago 
(ref. 10), then rediscovered by Stanton (ref. ll) for pipe flow, and later 
observed by Fritsch (ref. 12) for channel flow. For the special case of 
boundary- layer flow of an incompressible fluid along a flat plate, the 
analytical exploitation of these two similarity laws to predict skin 
friction and other important turbulent boundary- layer characteristics 
was initially given by Karman (ref. 13)* The fundamental basis of this 
velocity-profile similarity method, which is totally different in prin- 
ciple from the mixing-length method of Karman, has been discussed by 
Millikan in reference l4. A detailed review of the combined wall law and 
velocity-decrement law, including utilization of the equations of motion 
to calculate the shear-stress distribution across the boundary layer, may 
be found in the recent thesis of Coles (ref. 15) • 
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An active interest of aerodynamic is ts in the basic knowledge of 
turbulent boundary layers has been revived because of the necessity of 
extending methods of predicting incompressible boundary-layer character- 
istics to high Mach numbers where large variations in density occur across 
the boundary layer. Similarity laws for the variation of mean properties 
across a turbulent boundary layer in compressible flow have not as yet 
been found, although Coles has formulated a generalization of the wall 
law which involves an unknown function of Mach number to be determined 
from experiments. In the absence of detailed experiments to provide a 
sound foundation for analysis, recourse to mixing-length calculations and 
similar unproven hypotheses has been made in the past by a number of 
investigators in an attempt to obtain information about the turbulent 
boundary layer at high Mach numbers. The results, unfortunately, have 
been almost as varied as the number of investigators. This is indicated 
by figure 1 which shows the results (see refs. l6 to 31 ) of various theo- 
retical analyses of skin friction on an insulated plate in compressible 
turbulent flow.^ Since a brief discussion of these analyses was provided 
in a summary report of the present research (ref. 36 ) and also may be 
found in more detail in Coles^ thesis (ref. 15), no further discussion 
is given here. It will suffice, instead, to emphasize only what is clearly 
indicated by figure 1; namely, the present theoretical knowledge of tur- 
bulent boundary layers in compressible flow is severely limited. 

The practical reasons for requiring detailed knowledge of turbulent 
skin friction at high Mach numbers are twofold. The first concerns the 
obvious direct application to drag and, hence, performance estimations. 

The second, which is becoming of increasing practical importance, con- 
cerns the application to heat-transfer calculations, or the so-called 
"aerodynamic heating" problem. Although various estimates of heat- 
transfer rate (see Rubesin’s compilation in ref. 37 or the compilation 
of curves in ref. 36 ) are about equally as varied as those of skin 
friction, they, nevertheless, have one characteristic in common: a cal- 

culation of heat transfer first requires a knowledge of skin friction. 
Average skin-friction measurements can be used directly to estimate 
average heat-transfer rates for conditions of constant surface temperature. 
Such measurements also can be used indirectly as a guide in estimating 
local skin friction and local heat transfer, the two fundamental quantities 
used in practical calculations for conditions of variable surface tem- 
perature and pressure. 

In cognizance of the limitations of present turbulence theory and the 
necessity of knowing skin friction before either drag or heat transfer 

^Not indicated on this plot are the results of Dorodnitsyn (ref. 32 ) and 
and Kalikhman (ref. 33) which coincide with Karma’s estimation, and 
the result of Tetervin (ref. 3^) which predicts no effect of compres- 
sibility. Also not shown are the "power-law generalization" of Cope 
(ref. 18 ) which is tedious to compute at high Mach numbers because of 
poor series convergence, and the analysis of Donaldson (ref. 35) which 
apparently was not carried far enough to predict the effect of Mach 
number for constant- length Reynolds number. 
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can "be estimated, numerous measurements of turbulent skin friction in 
compressible flow have been made during the last few years. In this 
connection, reference is made to the measurements of Wilson (ref. 2l); 
Rubesin, Maydew, and Varga (ref. 25); Liepmann and Dhawan (ref. 28); 

Coles (refs. 15 and 39); Brinich and Diaconis (ref. 4o); Bradfield, 
DeCoursin, and Blvuner (ref. i<-l); Cope (ref. i;2); Weiler and Hartwig 
(ref. 43 ); Spivack (ref. 44); Monaghan and Johnson (ref. 45); Ladenburg 
and Bershader (ref. 46); Bloom (ref. 47); and Hakkinen (ref. 48). Results 
of these measurements are discussed later, especially those of Coles. 


The purpose of the present investigation was to measure the average 
friction for turbulent flow along the external surface of relatively long 
cylindrical bodies of revolution at both supersonic and subsonic speeds. 
For long cylinders the boundary-layer thickness becomes comparable to the 
cylinder radius, under which conditions appreciable departiires from flat- 
plate skin-friction values would be expected, according to the analyses 
of Jakob and Dow (ref. 49) , Landweber (ref. 50) , and Eckert (ref. 51)* 
Consequently, in the present tests various cylinders were investigated 
ranging in length from 8 to 23 diameters. The inclusion of subsonic 
measurements as an. integral part of the investigation was believed neces- 
sary for an accurate evaluation of the effect of Mach munber, since ade- 
quate skin-friction data were not available for incompressible flow along 
cylinders. It was considered desirable, further, to determine the fric- 
tion by measuring forces directly, rather than to deduce the friction 
indirectly from boundary-layer surveys according to the momentum method. 

As will be seen later, these two features proved quite valuable in mini- 
mizing several uncertainties common to most previous investigations . 

The present experiments were initiated in 1949, though the technique 
of measurement finally employed was not developed satisfactorily until a 
year later. The data presented herein were obtained during 1950-52. 

Some of the principal results have been simmarized previously in reference 

36. 
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NOTATION 


average skin-friction coefficient. 


force 

(l/2)p„Uoo^ (wetted area) 


local skin-friction coefficient. 


2Ty 

Poo^oo^ 


cylinder diameter 

friction force on cylinder of length I 
friction force on cylinder of length AZ 


Z 


length of cylinder over which friction force was measured 
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length of cylinder necessary to produce boundary layer of 
momentum thickness Qq at Reynolds number^ 

Vo 

corrected length of cylinder^ Z-fAZ 
Mach number 
static pressure 

base pressure (average of readings from five orifices) 
total pressure behind normal shock wave 
Reynolds 'number 

Reynolds number based on momentum thickness 

temperature 

velocity 

distance from surface 
boundary-layer thickness 

boundary- layer momentum thickness for cylinder. 



coefficient of viscosity 

N 

kinematic viscosity, 
mass density 
shear stress at wall 


Subscripts 


incompressible flow 

conditions at station 0 
(See fig. 2 and table I.) 

local conditions at outer edge of boundary layer 
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00 average of conditions at outer edge of boundary layer over 

length L 

w conditions at wall 

t total conditions 

APPARATUS AND TEST METHODS 
Wind Tunnels 


Tests were conducted in the Ames 1- by 3-foot supersonic wind tunnels 
No. 1 and No. 2. Wind tunnel No. 1 is of the closed-circuit, continuous- 
operation type equipped with a flexible-plate nozzle that provides a vari- 
ation of Mach number from 1.2 to 2.2. Total pressure in tunnel No. 1 can 
be varied to provide Reynolds numbers from 1 million to 10 million per 
foot. Wind tunnel No. 2 is of the nonretiirn, intermittent-operation type 
also equipped with a flexible -plate nozzle to provide a variation of 
Mach number from 0.5 to 3.8. Compressed air at a majcimum pressxire of 
6 atmospheres is obtained from a reservoir and is expanded through the 
nozzle to atmospheric pressure. The total pressure is controlled by means 
of a throttling butterfly valve between the reservoir and the tunnel 
settling chamber, providing a variation of Reynolds number from 7 million 
to 20 million per foot. 

Primarily because of the throttling method employed, the tunnel tur- 
bulence level in the No. 2 tunnel is much higher than that in the No. 1 
tunnel. Some unpublished hot-wire measurements made by Howard A. Stine 
of the NACA Ames Laboratory have indicated the turbulence level in the 
No. 1 tunnel to be between 2 and 3 percent as measured in the reservoir 
at a location where the velocity was of the order of 15 feet per second. 

In the one-year interval between the last of the shin-friction measure- 
ments and these recent hot-wire measurements, the No. 2 tunnel was com- 
pletely revised. Consequently, recent hot-wire measurements in the No. 2 
tunnel, which indicate the turbulence level at the time of this writing 
to be between 2 and 10 percent, depending on operating conditions of 
total pressure and Mach number, do not apply to data presented in this 
report. As a result of the revisions, though, it is known that the 
transition Reynolds number on a 10*^ cone (measured by the method described 
in ref. 52) was increased several fold. It is believed, therefore, that 
the turbulence level for the No. 2 tunnel as it existed during the skin- 
friction measurements was considerably greater than that for the revised 
tunnel. It will be shown later that the difference in turbulence levels 
of the two wind tunnels had no measirreable effect on turbulent skin fric- 
tion at a Mach number of 2 . 
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The water content of the air in both wind tunnels was maintained at 
less than 0.0003 pound of water per pound of dry air. Consequently, no 
correction was made for humidity effects. 


Force Measurements 


The friction force f acting along a cylinder of length I and 
diameter D was determined from two drag measurements by subtracting 
the absolute foredrag (sum of measured drag force and the product 
Pb^D^/4) of a short cone-cylinder of length Z» from that of a much 
longer cone-cylinder of length Z+Z • . (See fig. 2.) Photographs of two 
models employed are shown in figure 3. Models were anodize-f inished 
aluminum and appear black in the photographs as contrasted with the 
cadmium-plated shrouds. For each model the cylinder diameter was 1 inch 
and the included cone angle was 20°. Three different lengths Z+Z' were 
investigated: 10, 15, and 25 inches, as measured from the tip of cone 

to base of cylinder. In the final Cp data to be discussed below, these 
lengths corresponded to cylinder fineness ratios of approximately 8, 13, 
and 23, respectively. It was found necessary to provide boundary-layer 
trips in order to assure a turbulent boundary layer in all but a few 
test conditions. The trips employed will be discussed later. The magni- 
tude of the friction drag was between 20 percent and 50 percent of the 
over-all body drag and is believed to have been measured with a correspond- 
ing accuracy of from U percent to 2 percent, depending upon length- 
diameter ratio and test conditions . 

Total drag was measured by means of a three -component electrical - 
strain-gage balance which had been -specially adapted to the measurement 
of axial force for the present investigation. The balance and supporting 
stings were shielded from the oncoming flow by shrouds in the manner 
shown by the photograph of figure 3. The cylinders were alined with the 
shrouds by the use of ball-bearing guides through which the stings were 
passed. 

Base drag was measured by means of five base-pressure orifices 
connected to a multimanometer using dibutal phthalate as the manometer 
fluid. Three of these orifices were located evenly spaced in the face 
of each shroud, the fourth being located on the sting, and the fifth 
inside the balance chamber. In addition to determining the base drag, 
the five pressure measurements provided a sensitive indication of aline- 
ment of cone -cylinders with shrouds and also confirmed the absence of 
leaks in the balance -shroud enclosure. No data were taken unless the 
five base pressures were in essential agreement. A reference tube of the 
manometer was connected to a vacmmi pump during supersonic testing but 
was left open to the atmosphere when operating at subsonic speeds . A 
McLeod gage was used to measure the pressure of the vacuum reference. 
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Apexes of the two cone-cylinders needed for one friction-force 
measurement were made to occupy the same position in the tunnel test 
section, thereby eliminating the need of an axial buoyancy correction. 

This technic[ue also enabled axial pressure distribution along the cylinder 
over which the friction force acted to be measiared by eight equally spaced 
orifices located in that portion of the shroud surface which corresponded 
to the cylinder of length i. Static-pressure distributions obtained in 
this manner were foiind to agree with those measured on a similar cone- 
cylinder model without a gap. 

Allowance was made for the fact that the two drag measurements used 
in determining the friction force f were made during two separate runs 
under slightly different conditions of tiinnel pressure and temperature. 
Faired curves of the absolute foredrag coefficient of the short cone- 
cylinder plotted against Reynolds number per unit length were used to 
determine the absolute foredrag at the particular Reynolds number cor- 
responding to the drag measurement of the long cone-cylinder. 

The drag gage was calibrated prior to each day' s testing, and a check 
of the balance was made before and after each test to assure that the 
friction in the balance system was negligible while measuring forces. 

This practice showed that periodic cleaning of the balance was essential 
to the satisfactory measurement of the friction force. Further uncertain- 
ties were minimized by repeating each set of measurements at least once. 


Method of Determining Effective Reynolds Number 
and Corrected Skin-Friction Coefficient 


As indicated in the center sketch of figixre 2, the boundary- layer 
thickness at the leading edge of the cylinder over which the friction 
force f was measured was not zero but some value Sq- Hence, to provide 
a common basis of comparison, it was necessary to apply corrections in 
order- to convert the force measurements to conditions of completely tur- 
bulent flow starting with zero boundary -layer thickness at the leading 
edge of a cylinder. For this purpose one can imagine the cylindrical 
surface to be extended upstream a short distance Al (see bottom sketch 
in fig. 2) such that the turbulent boundary- layer momentum thickness 
developed over the length Al just matches the momentum thickness 9q 
corresponding to the boundary- layer thickness 6 q. Boundary-layer surveys 
in the plane of station 0, as illustrated in the sketch second from the 
bottom in figure 2, were used to determine boundary- layer profiles from 
which Al and the corresponding friction force Af were calculated. 

The surveys showed that Qq depended on Moo^ Ro? nose shape, etc., in a 
complicated manner that could be determined only by measurement. As a 
result, it was necessary to obtain boundary -layer surveys in this position 
for all the conditions under which the direct-force measurements were 
obtained. This proved to be a time-consuming task since it involved 
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making a total of about 270 boundary-layer surveys. The small friction 
force Af that would act over the length AZ was readily determined 
from the equation 


Af = JtDPoUo^Qo 


which represents conservation of momentum. The length AZ was detennined 
from an equivalent form of the momentum equation 


AZ 


20 


o 




Actually, the skin-friction coefficient Cp , based on the Reynolds number 

UqAZ/Vo, was not known until AZ was known. Hence, the method of deter- 
mining and AZ was, in principle, an iterative one, using the 

equation 

CFo ■ CFio (^) 


with evaluated at a Reynolds number UqAZ/vq from the Karman- 

Schoenherr equation 

0.242 


Zogio (ReCp^) 


in conjunction with the ratio (Cp/Cp^) evaluated at an arbitrary Reynolds 
number from a plot of uncorrected force data (this latter ratio does not 
depend significantly on Reynolds number) . It turned out that the cor- 
rections AZ and Af were sufficiently small and the correction method 
so rapidly convergent, that, after the first few trials, iteration was 
found to be unnecessary. 

The average skin-friction coefficient Cp and corresponding Reynolds 
number Re employed throughout are based on the total length L = Z+AZ 
and the total force f+Af according to the defining equations 

f + Af 

(1/2) p„u„ 2 [nD(Z + AZ)] 

U»( Z + AZ ) 


Reference quantities such as sind correspond to conditions 

at the outer edge of the boundary layer averaged over the length L. These 
were determined from measurements of total pressure, static pressure. 


10 


NACA TN 3097 


total temperature, and the Sutherland equation for viscosity (with 
S = 216 and = O. 36 I x 10“® Ih sec ft”^ at 500^ R), assuming the air 
to he a perfect gas. 

Since it is not possible in any experiment to obtain fully developed 
turbulent flow starting precisely at the leading edge, some assumption 
inevitably is involved in any method of determining the effective starting 
position of turbulence, even though the method appears to be based on 
measured quantities such as momentum thickness. The method described above 
is based on the assumption that the ratio does not depend sig- 

nificantly on Reynolds number. This assumption should be reasonable for 
the Reynolds numbers of the investigation. (The lowest Reynolds number 
at which the starting-length correction was applied was 0.4 x 10®.) Other 
methods have been employed which involve different assumptions (see refs. 
15 , 21, 25 ^ 45 , and 46). None of the methods is exact. Hence, in order 
to obtain accurate and consistent data, it is necessary that the starting- 
length correction be sufficiently small so that various correction methods 
would yield substantially the same end result. Such conditions were met 
in the present experiments inasmuch as the starting-length correction to 
skin-friction coefficient was relatively small. This correction generally 
decreased with increasing Mach number, varying between 10 and 3 percent 
for l/d =8, 6 and 2 percent for l/d = 13^ and 4 and 1 percent for 
l/d = 23* In the most extreme case, a 10 -percent error in the correction 
itself would introduce 1-percent error in Cp. 


Boundary-Layer Survey Apparatus and 
Reduction of Survey Data 


The apparatus used in obtaining all survey measurements is shown in 
the photograph of figure 4 which illustrates the position of the various 
components of the equipment while making boundary- layer surveys in the 
plane of station 0. Dimensions locating this station, as well as the 
other stations at which surveys were made, are given in table I. The sur- 
vey instrument consisted essentially of a housing pinned to a supporting 
member. A boundary- layer probe which was attached to the housing could 
be moved in a vertical direction since the housing could be rotated 
approximately 2^ in the vertical plane passing through the tunnel axis. 
Because the probe opening was located about 15 inches upstream of the 
pivot pin, the probe mouth could be considered to move in a straight line. 
Motion was provided by the use of an electric motor located outside the 
wind tunnel and connected to the survey instrument by a flexible shaft. 
Speed reduction was such that position in the vertical direction could be 
controlled to about 0.0005 inch. 

Two probes were used, the dimensions of which are given in the table 
of figure 5- These probes were constructed of stainless-steel hypodermic 
tubing following the spirit of the suggestions of reference 53- Time lags 
of the order of 5 seconds for the highest Mach number and a fraction of 
a second for the lowest Mach mamber were obtained. 
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Probe -opening position and pitot pressure were determined by an 
electrical-resistance strain gage and a pressure cell, respectively, 
both of which were located within the instrument housing* The base of 
the strain gage was fixed to the supporting member in such a manner that 
deflection of the gage was increased by downward movement of the housing* 
Characteristics of the pressure cell employed can be found in reference 
5^* Static pressure was determined (with the probe outside the boundaiy 
layer) by means of orifices located in the surface of the survey model at 
the various boundary-layer survey stations. At each station two orifices 
were located in a vertical plane on opposite sides of the model from each 
other in order to indicate the degree to which the model had been alined 
with the stream. This precaution, however, did not assure obtaining iden- 
tical boundary-layer sui^eys from the top and bottom of the survey model 
except in a few cases due to the high sensitivity of boundary-layer pro- 
file to small stream angles. 

Two methods were used to determine the point at which the probe con- 
tacted the surface of the survey model in wind tunnel No. 1* Since the 
survey model was electrically insulated from the survey instrument, the 
point at which contact occurred should have been indicated by a marked 
decrease in the resistance between probe and model as measured by an elec- 
trical analyzer. Actually, either the probe or the model vibrated some- 
what so that approaching contact was at first indicated by intemittent 
indications of the meter and, as the probe continued toward the model, 
the mean value of the indicated resistance steadily decreased until solid 
contact was finally made* In order to determine the value of the resist- 
ance which corresponded to effective probe contact, the region near con- 
tact was viewed from outside the wind tunnel through the telescope of a 
cathetometer . The combined methods enabled the contact position to be 
determined with an estimated accuracy of about 0*001 inch for measurements 
in the No. 1 wind tunnel. In the No. 2 wind tunnel, the optical method 
could not be used and only the electrical method was employed for the 
contact-position determination. Consequently, the estimated accuracy of 
this measurement in the No. 2 wind tunnel was of the order of 0.003 inch. 
Inasmuch as the contact position was used as the zero point of the bottom 
surface of the probe, the absolute position of the vertical scale is uncer- 
tain by this amount. At worst, such an uncertainty could mean a 20-percent 
error in 0 q, with a corresponding error of 2 percent in Cjp. The effect 
of probe or model vibration on the measurement of pitot pressure is not 
known. 

Computations necessary to obtain boundary-layer velocity profiles 
and other boundary-layer characteristics from the data obtained during 
surveys were carried out on a card-programmed electronic calculator with 
the usual assumptions of Prandtl number equal to unity and constant 
static pressure through the boundary layer. One difficulty ordinarily 
encountered with boundary-layer computations by these machines is that 
a long table search is required for obtaining Mach number from the ratio 
of static pressure (p) to pitot pressure (p^O* This step was simplified 
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somevhat by employing 23 straight-line segments as an approximation to 

the curve of (p^*-p)/|pu^ as a function of p/p^' • Coordinates of the 
end points of the straight-line segments are given in table II, which 
covers the range of Mach numbers from zero to infinity. 


RESULTS OF MEASUREMENTS TESTING VALIDITY 
OF EXPERIMENTAL METHOD 


Early in the course of measurements it became apparent that the 
uncorrected skin-friction data depended significantly on such things as 
the size and location of boundary-layer trip, the length of nose ahead 
of the cylindrical measuring section, and the particular wind tunnel in 
which the data were obtained. This situation arose primarily because the 
location of transition - and, hence, the effective origin of fully devel- 
oped tiirbulence - varied with changes in each of these conditions. It 
was hoped that if the effective origin of turbiilence were determined for 
each test condition, and if the production of completely turbulent flow 
by a boundary-layer trip were assured, then the final corrected data would 
not be sensitive to the techniques employed in obtaining them. Numerous 
supplementary tests, which are discussed below, were made in order to 
investigate such matters and thereby to determine the degree to which the 
final data depend on the particular experimental methods employed. 


Tripping Devices to Insure Turbulent Flow 


Surface-probe pitot-pressure surveys with various trips applied were 
made to find a satisfactory method of obtaining completely turbulent flow. 
Dimensions of the various trips investigated are given in figure 6. The 
results for M„ = 2.9 and for the lowest tunnel pressure employed in the 
No. 2 tunnel at this particular Mach n\imber are shown in figure 7 . The 
Mach number (m) ^.^ ference ^sed in the ordinate of this figure is the local 
value that would exist in inviscid flow, as determined by the method of 
characteristics. It may be seen that transition on the smooth body appar- 
ently occurred between about 4 and 7 inches from the cone tip. Conse- 
quently, some device, such as a boundary-layer trip, wa.s needed to obtain 
a completely turb\ilent boundary layer over the cylinder. As figure 7 
shows, trip 6 moved transition upstream some but not enough. Trip 5 
appears sufficient, and trip 1, which is located the farthest upstream, 
appears more than sufficient \inder these particular conditions. 

Although trip 5 was adequate at M„ = 2.9, it was not adequate at 
Moo = 3*6. Trip 1 proved to be adequate for every Mach niimber and Reynolds 
number investigated in the No. 2 wind tunnel; hence, it was selected for 
use in all final measurements of turbulent skin friction conducted in that 
particular tunnel. In the No. 1 wind tunnel, however, the turbulence level 
and Reynolds niimbers obtainable are much lower, and it was found necessary 
to employ two trips in tandem (trips 2 and L proved adequate) . 
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As an incidental note of possible interest, the point of beginning 
fuzziness in the boundary layer on various schlieren photographs (not 
shown) from the No. 2 tunnel was in all cases close to the apparent end 
of transition indicated by surface-pressure surveys. This correlation 
is believed not to be a general one, but dependent on the particular 
wind tunnel, since fair correlation was observed in the No. 1 tunnel with 
the apparent beginning of transition indicated by surface -pres sure surveys. 

An independent demonstration of the state of boundary-layer flow at 
the beginning station of measurement (station O) was provided by the 
boundary-layer velocity-profile sui-veys. Typical examples from the No. 1 
tunnel are shown in figure 8. With the addition of a trip, the actual 
velocity profiles, (using y as ordinate, fig. 8(a)) were thicker than for 
the smooth body; whereas the dimensionless velocity profiles (using y/S 
as ordinate, fig. 8(b)) were much fuller than corresponding profiles 
without a trip. The fuller profile in (u/ue, y/e) coordinates is more 
characteristic of turbulent profiles. Once the boundary layer was ade- 
quately tripped, further additions of tripping devices only increased 
the boundary-layer thickness, but did not alter its dimensionless velocity 
profile. Typical data from the No. 2 tunnel which demonstrate this latter 
point are shown in figure 9(Q') ^ and the results of measurements of skin- 
friction coefficient with various trips and without a trip are shown in 
figure 9(b) . A curve is not faired through the transition-region data 
for the smooth cone-cylinder (circle points in figure 9(b)), inasmuch as 
these data are believed to exhibit primarily such things as the effects 
of turbulence level of the No. 2 tunnel. Also, the particular method 
employed to determine the effective boundary-layer origin is not applica- 
ble to transitional or laminar flow. It is apparent, though, that the 
differences in skin-friction coefficient between various trips are not 
large, and that the two data points for natural transition at the highest 
Reynolds numbers (near 20 million for l/d = 13) agreed well with the 
data obtained for tripped boundary layers. 

A comparison of the velocity profile for the smooth body in the 
No. 1 tunnel (fig. 8(b)) with that for the same body in the No. 2 tunnel 
(fig. 9(a)) shows marked differences. In both cases the Reynolds number 
based on distance from cone tip to station 0 is about 2 million. Since 
transition on a smooth 10 ^ cone at M« = 2 begins at a Reynolds number 
of 3»1 million in the No. 1 tunnel, and at 0.7 million in the No. 2 tunnel 
(see ref. 52), it is to be expected that a laminar-type profile would be 
observed in the No. 1 tunnel and a transitional or near -turbulent -type 
profile would be observed in the No. 2 tunnel. Such expectations are 
borne out by the data of figures 8(b) and 9(a). 

A point that should be noted from figure 8 is the appreciable differ- 
ence sometimes obtained between velocity profiles at the top and bottom 
surfaces of the cylinder. This difference is due to a very small mis- 
alinement (a few tenths of a degree) of the model axis with the undisturbed 
stream direction. The high sensitivity of velocity profiles to small 
angles of attack on bodies of revolution has been noted before by several 
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investigators^ and indicates that the circumferential variation in 0 
should he determined for such bodies when measuring skin friction by the 
momentum method. Such variations do not affect determinations of Cjr 
by the direct-force method employed in this investigation because measure- 
ments of 0 affect only the small correction to the force data. More- 
over^ such variations do not directly affect the skin-friction force 
measurements^ since these were found to be independent of angle -of -attack 
variations up to the order of ± 1 ^. 


Effect of Nose Shape 


Two different nose shapes were employed at Mqo = O.8I to determine 
whether or not the subsonic data depended materially on this parameter. 
Although the measured momentum thickness at station 0 differed by a 
factor of approximately 2 for the two nose shapes shown in figure 10 , the 
final skin-friction data were essentially independent of nose shape. 

(See fig. 10 .) This should not be surprising in view of the fact that 
the only appreciable alterations in pressure distribution brought about 
by varying nose shape are confined to a small region downstream of sta- 
tion 0 ; hence, the main effect on average skin-friction coefficient is 
believed due to the differences in 0q, which can be corrected for, just 
as for the case previously mentioned of different boundary-layer trips. 


Attainment of Equilibri-um in Blowdown Tunnel 


Since tiinnel No. 2 is of the blowdown type, the measurements in 
that tunnel were always made under conditions of slight heat transfer 
to the air stream. Figure 11 (a) shows a typical time history of the 
wall temperature T-^ expressed as a fraction of the reservoir total 
temperature Tt (which decreases at the rate of about F per minute). 
Clearly, after about 5 minutes Tw/^t became constant. The limiting 
value, 0.940 at = 3 . 4 , implies a recovery ratio (Tv-Too)/(Tt "'^00) 
of approximately O.9I which is slightly higher than the usual recovei^ 
factor (0.8910.01) because of the small amount of heat transfer present 
even after Tv/T-^ became constant. This small heat transfer (wall tem- 
perature within 2 percent of recovery temperature) is believed to have 
a negligible effect on skin friction. Only one set of measurements was 
made of Cjr as a function of time, with the results shown in fig- 
ure 11(b). The small variation probably is not significant, since 
2-percent variation is within the limit of experimental error. 
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Comparison of Data From Wind Tunnels 


At Mach numbers near 2.0, measurements could be made in tunnels 
No. 1 and No. 2 at overlapping Reynolds numbers. This provided an oppor- 
tunity to compare data talien under considerably different conditions. 

As noted earlier, the No. 1 tunnel is a relatively lov-turbulence -level 
tunnel (Re for transition on 10^ cone is 3.1 x 10® at = 2.0) of 
continuous -operation type in vhich the models were in thermal equilib- 
rium. The No. 2 tunnel is a relatively high-turbulence-level tunnel 
(Re for transition on lOO cone is 0.7 X 10® at = 2.0) of the blow- 
down type in which the models were releasing small amounts of heat to 
the stream, as discussed above. These combined differences, however, 
do not affect the observed velocity profiles at = 2, as is evident 
from figure 12(a). ^ Likewise, these differences have no appreciable 
effect on the measured skin-friction coefficient at = 2, as is evi- 
dent from figure 12(b). It is not known whether the subsonic measure- 
ments, obtained only in the No. 2 tunnel, are affected by the relatively 
high turbulence level. 


Effect of Pressure Gradient 


At a Mach niimber of 2.0, the flexible -plate nozzle walls pui*posely 
were distorted in order to obtain friction data with various pressure 
distributions along the cylinder. Two such distributions, from a total 
of three, are shown in figure 13(a) plotted in the form of Mach number 
distribution. The corresponding skin-friction data are shown in fig- 
ure 13(b). Evidently, the average turbulent skin-friction coefficient 
at this Mach number is not significantly affected by variations in pres- 
sure of the amounts investigated. It may be noted that the third dis- 
tribution investigated (not shown) yielded the same skin friction to 
within about 2 percent as the two distributions shown. It is possible 
that the effect on local skin friction at individual points along the 
cylinder, which could not be measured, was considerably greater than 
the small integrated effect on average skin friction. It is also pos- 
sible that the observed insensitivity to pressure gradient does not 
persist to Mach numbers far removed from 2.0 because the coefficient of 
the pres sure -gradient term in the momentum equation for boimdary layers 
happens to vanish for turbulent flow at Mach numbers of about 2.4 (see 
ref. 55 or I 5 ) . 

Mention should be made of the fact that the two pressure distribu- 
tions shown in figure 13 (a) differ from each other by an amo-unt greater 
than the difference between the normal pressure distributions of the 

^Both velocity profiles shown in figure 12 (a) were measured at the top 
surface of the survey model. In each of these cases, however, profiles 
obtained at the bottom surface agreed with those shown. 


l6 


NACA TN 3097 


present investigation and a constant pressure. This difference in 
pressure distribution is comparable to the differences between pressure 
distributions for the various experiments on flat plates to which com- 
parison with the present results will be made subsequently. Also, it 
is noted that the two pressure distributions shown differ from eacb 
other by roughly the same amount that the pressure distribution on a 
3-percent-thick biconvex airfoil differs from a constant pressure; or, 
alternatively, by roughly the same amount that the pressure distribution 
on a parabolic -arc body of revolution of fineness ratio 15 differs from 
a constant pressure. The range of pressure variation covered by the 
present tests, therefore, is of some practical significance. Since no 
important effect on Cp of variations in pressure distribution was 
observed, the average skin-friction data to be presented later are con- 
sidered applicable not only to the ideal case of constant pressure but 
also to thin pointed airfoils. Such data, though, apply only approxi- 
mately to bodies of revolution not having a constant diameter since, in 
such cases, corrections should be made to allow for the relative thin- 
ning of the boundary layer over an expanding nose (Mangier type, correc- 
tion) and for the departure from two-dimensionality over an aftersection 
of the body due to the boundary -layer thickness becoming comparable to 
the body radius (l/D-type correction). 

Although the departures from uniform pressure distribution did not 
affect average skin friction substantially, they did have noticeable 
effect on the velocity profiles over a small portion of the cylinder 
downstream of the cone -cylinder juncture. This is illustrated in fig- 
ure Ik which shows five profiles at various stations along the cylinder. 
In the region near the shoulder (station O) where the pressure varies 
most rapidly, the velocity profile is fuller than at downstream stations. 
This effect does not persist very far downstream, as may be seen from 
the figure. Little can be said about changes in dimensionless velocity 
profiles downstream of station 2 except that they were not greater than 
the random differences often observed between top and bottom surveys at 
a given station. Similarly, the rate of growth of boundary-layer thick- 
ness along the cylinder cannot be accurately specified because of this 
sensitive effect of model misalinement on the velocity profiles. 


FINAL RESULTS AND DISCUSSION 

Skin Friction for Various Values 
of l/D, M„, and Re 


The principal results of the investigation showing Cp as a func- 
tion of various parameters are presented in figures 15^ I6, and 17. 

These figures will be discussed together in view of the interdependence 
of the various data. The fairing of curves in figure 15 has been influ- 
enced by the data of figure I6. In figure 16 the ratio Cp/Cp^ for a 
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given Reynolds number is plotted as a fiinction of Re. The data shown 
are typical of all data obtained. In some cases, the data points indi- 
cate a small decrease in Cp/Cp^ with increasing Re (e.g., fig. l6(d)); 
in other cases, a small increase would be indicated (e.g., the data for 
l/d = 13 at M«o =2.5 in fig. 15 would indicate such a trend if Cp/Cp. 
were plotted as a fimction of Re); but in most cases, no systematic 
variation in Cp/Cp^ with Re is indicated. The variations sometimes 
observed are small enough to be within experimental error. Consequently, 
the curves in figure 15 have been faired to be consistent with the result 
of figure 16, indicating no substantial effect of Re on Cp/Cp^ over the 
Re range covered by a given model. If the Re range were increased 
many times, a small effect on Cp/Cp^ might be expected. 

In figure 15, as well as in all other figures where it is not speci- 
fied otherwise, the reference skin-friction coefficient Cpj^ for zero 
Mach number arbitrarily is defined as the value given by the Karman- 
Schoenherr equation for flat -plate flow (Cp^) . This arbitrary selec- 

tion from the many semi empirical equations for Cp^ was made primarily 
because, at the Reynolds numbers for which measurements in compressible 
flow are available, it provides results intermediate to the res-ults of 
other careful evaluations of incompressible turb\ilent skin friction. 

(See Appendix for comment on this point.) Discussion is presented later 
as to a more proper reference Cp^ for the present experiments on 
cylinders . 


Plots of Cp/(Cp^) as a function of are presented in figure 

17(a) and 17(b) for the^tlree l/d values investigated. As noted in the 
figiire subtitles, figure 17(a) presents individual data points for certain 
nominal Reynolds numbers; whereas figure 17(b) presents values picked 
from curves faired through measurements at all Reynolds numbers investi- 
gated. Both figures show no large effect of cylinder length on C-p. 

Since the maximum effect observed is only about twice the estimated experi- 
mental error, it is not considered significant that a few data points, 
especially in figure 17(a), appear inconsistent with the main body of 
data (fig. 17(b)) which exhibits the trend of increasing skin friction 
with increasing length -diameter ratio. The increment in Cp between 
l/d = 8 and l/d = 23 appears roughly consteint for most Mach numbers, 
suggesting that perhaps at Mach numbers of the order of 5 or higher a 
substantial percentage effect of l/d on Cp might exist. This trend 
would be consistent with the increase in boundary -layer thickness at the 
higher Mach numbers, but it is by no means certain, in view of the fact 
that the differences observed are only about twice the estimated experi- 
mental error. 


With regard to theoretical calculations of the effect of l/D on Cp, 
it is to be noted that the analyses of Jakob and Dow (ref. 49), Landweber 
(ref. 50), and Eckert (ref. 51) predict an increase in Cp with increas- 
ing l/d. The analysis of Jakob and Dow, which assumes that the rate of 
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"boundary -layer grovth along a cylinder is identical to that along a 
plate ^ predicts approximately 15-percent increase in friction for 
l/d = 23 as compared to l/d = 8. The analyses of Landweher and Eckert, 
which assumed that the relation between friction coefficient and boundary- 
layer thickness is identical for cylinder and plate, predict only 1.5“ 
percent increase. The corresponding experimental result from the sub- 
sonic measurements is approximately a 4-percent increase. Unfortunately, 
the basic assumptions of these two analyses cannot be evaluated from the 
present experiments because of the previously mentioned difficulty of a 
circumferential variation of boundary-layer thickness often encountered 
at a fixed longitudinal station. 


Thus far, the reference value employed for of cylinders is the 

Karman-Schoenherr value for flat-plate flow. It is now apparent that the 
correct reference value for Cp. will depend somewhat on the l/d ratio 
of the cylinder. Consequently, the true effect of Mqq on Cp for a given 
l/U would have to be evaluated by using values of Cp^ appropriate for 
each value of l/D. Sufficient subsonic measurements were made to provide 
such an evaluation. An experimental value (Cp^^) of skin-friction 


exp 


coefficient for each cylinder at = 0 was determined from the two sets 

of subsonic measurements by applying a correction to allow for the small 

effect of compressibility at subsonic Mach members. The correction was 

determined by averaging the predicted values of Cp/Cp^ given by the 

equations of Cope, Tucker II, and Wilson, and amounted to a 2.5- and 6.1- 

percent increase in the data for Moo =0.51 and Moo = 0.8l, respectively. 

The two values so obtained were averaged for each cylinder. Values of 

(Cp^O deduced in this manner from the data of figure 17(a) for l/d 
^ 1 exp 

ratios of 8, 13 , and 23 are 2.5, 0.2, and ^-0.6 percent lover, respec- 
tively, than the value predicted by the Karman-Schoenherr equation. The 
corresponding corrections from the data of figure 17(h) are 2.2, 0.2, 

and -l.i<- percent, respectively. With (Cp . ) so determined, plots of 

1 exp 

Cp/CCp.) were made. Results from the-- data of figure 17(h) are shown 

in figure 17(c). There is little scatter in this latter plot (mean devi- 
ation from faired curves is about 1 percent) and no clearly discernible 
effect of l/d. Consequently, data points for the various fineness 
ratios were averaged at each Mach number. The resulting values, tabu- 
lated in tables III (a) and Ill(b), corresponding, respectively, to data 
of figures 17(a) and 17 (l), can be compared directly with flat-plate 
data, in view of the absence of any appreciable l/d effect on 

Cp/(Cp.) . The values in table Ill(b) are perhaps more representative 
1 exp . , N 

of the main body of experimental data than the values in table III(a;, 

although the differences are not significant. Table III (a) was published 
previously in reference 36 . 
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Comparison of Results With Other Investigations 

In figure I8 tvo velocity profiles from the present investigation on 
cylinders are compared with a velocity profile from Coles* measurements 
(ref. 15) on a flat plate. These profiles are for Moo = 2.0 and for 
Reynolds numbers near 6 million. It is seen that there is an appreciable 
ciillerence between the profiles for plate and cylinder^ the cylinder pro- 
file being the fuller of the two. 

A compilation of results from the several experiments in which shin 
friction was determined by direct-force measurements is presented in 
figure 19.^ Here^ again^ the ratio of compressible to incompressible 
skin friction is plotted as a function of Mach number. The Reynolds 
numbers of the three experiments are sufficiently comparable, in view of 
the results shown in figure 16, so that the indicated effect of Mach num- 
ber is not obscured by such differences in Reynolds number. As noted in 
figure 19^ the Karman-Schoenherr has been used as reference for the 

data of Coles and Hakkinen;"^ whereas experimental values of incompres- 
sible friction have been used as reference for the present data. Because 
of the use of a different reference c^^ and a different effective 
Reynolds niimber (see Appendix), the data points representing Coles' experi- 
ments in figure I9 are appreciably, though perhaps not significantly, 
different from points that would be obtained by plotting values tabulated 

^The direct-force measurements of Bradfield, DeCoursin, and Blumer 

(ref. 4 l) on cones are not included in this figure. It is noted, how- 
ever, that these data for cones indicate higher friction than the data 
in figure 19 * The direct-force data of Liepmann and Dhawan (ref. 38) 
and Weiler and Hartwig (ref. ^ 3 ) ^ Tor which the effective origin of 
turbulence was not determined, also are not shown in figure I9. 

In an earlier paper (ref. 36) Coles' uncorrected data of refer- 
ence 39 were used for comparison, inasmuch as his corrected data were 
not available at the time of that writing. In figure I9 Coles' data 
have been referred to an average effective Reynolds number for reasons 
^ explained in the Appendix. 

Although Hakkinen made subsonic measurements of local friction, he did 
not determine the effective Reynolds number for these low-speed data. 

His effective Re for the supersonic data is based on distance from 
the position of maximum shear, which he determined. 

The value of c^^ corresponding to the Karman-Schoenherr equa- 
tion is ^ 

0.557 Cp 

^fi = 7 = 

0.557 + 2yCp^ 

where Cjp^ is as defined previously. 
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in reference 15* Consequently^ without some comment the apparently 
excellent agreement of Coles ^ c^ data with the present Cp data might 
be misleading. For example. Coles* data points would be lowered about 
2 percent if the so-called Prandtl-Schlichting equation for Cf^ were 
used, but would be raised about 3 percent if the equation proposed by 
Coles for cf^ were used. Similarly, his data points would be either 
raised or lowered the order of k percent, depending on the values of 
effective Reynolds number assigned to his data. Thus, the unusually 
close agreement indicated in figure 19 can be taken as meaning for cer- 
tain only the order of ±5-percent agreement between Cf/cf^ and Cp/Cp^, 
although the true agreement may possibly be as close as the figure 
suggests. 

Beyond a Mach number of about 5^ the absence of experimental data 
renders the estimation of skip friction increasingly uncertain as the 
Mach number is progressively increased. Perhaps a measure of the mini- 
mum uncertainty involved would be that indicated by the difference 
between the estimates made by Cope and Tucker. Although it is not shown 
in figure I9, these two estimates do not differ greatly from the experi- 
ments below Moo = 5* They diverge considerably from each other above 
Moo = 5^ as illustrated in the figure. Cope and Tucker each have advanced 
two estimates; the particular curves shown in figure I9 represent Cope’s 
”log-law” analysis and Tucker’s "mean-reference -temperature” hypothesis^ 
both of which have been recomputed for the common basis of a recovery 
factor of 1, a 0.76-power variation of viscosity with temperature, and 
the Kaman-Schoenherr equation for incompressible flow. It may be noted 
that several other analyses (e.g., Wilson and Frankl-Voishel) also agree 
fairly well with the skin-friction measurements below = 5^ and that 
many other equations easily could be constructed which would do likewise# 
Consequently, without experimental confirmation of the basic assumptions 
involved, the agreement in end result of any one of the various analyses 
with experiment does not necessarily represent any more than woiild be 
represented by agreement of an empirical equation with experiment. 

In figure 20 a comparison is made between the mean curve from the 
direct-force data of figure I9 and numerous other measurements (the 
effective Reynolds numbers for which are listed in table V) wherein the 
average skin friction was determined indirectly from boundary-layer 
surveys using the conventional momentum method.^ The agreement is not 
good. Except for the present measurements, these various momentum- 
method data have all been referred to the Karman-Schoenherr equation 

^Experiments wherein the virtual origin of turbulence, and, hence, the 
effective -length Reynolds number, was not determined experimentally are 
excluded from figure 20. Measurements in this category are, for example, 
the tunnel-wall measurements of Cope (ref. 42) and Weiler and Hartwig 
(ref. 43). Also, the preliminary measurements of Bloom (ref. 47) have 
been excluded because of possible condensation effects in his experi- 
ments. 
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for Cp^. Excluding one point, the many momentum-method determinations 
indicate higher friction than the direct-force data. This particular 
point, from the present experiments at = 3.6, represents a single 
determination in the blowdown (No. 2) tvuinel at a single Reynolds number. 
The accuracy of this point is not known. In contrast, the point from 
the present experiments plotted at Moo = 2.0 represents the average of 
many determinations on two different models in the continuous -operation 
(No. 1) tunnel at various Reynolds niunbers between 3 million and 10 mil- 
lion, all determinations of which showed from 5- to 10 -percent -higher 
friction than the direct-force data. For these measurements on cylinders, 
an average of the surveys on top and bottom of the cylinder was employed 
in all cases, with the friction force deduced from the difference in 
surveys at stations 0 and 7* Similar s-urveys at intermediate stations 
were employed to evaluate the small correction for pressure gradient in 
the momentum integral equation. 

There are at least two additional reasons for possible errors in the 
conventional momentum method of determining skin friction, even for flows 
wherein the pressure is constant and transverse secondary flow does not 
exist. First, the reading of a pitot tube near the surface where the 
greatest momentum decrement occurs is not entirely certain because of 
probe interference and perhaps vibration; second, the conventional momen- 
tum method neglects the contribution of the Reynolds normal stress to the 
momentum integral. Estimates of this latter contribution made by Ross 
(ref. 56), for the case Ma, = 0, indicate that the conventional momentum 
method would yield values about 2 percent too high. It is noted that the 
correction factor of Ross involves principally the shape parameter (ratio 
of displacement to momentum thickness) which greatly increases as the 
Mach number increases. 

It should be noted that a given error in evaluating the effective 
Reynolds number for momentum-method Cp data results in much greater 
final error than for direct-force Cf data. This situation arises 
because Re enters directly into the ordinate for momentum-method deter- 
minations of average friction (Cp = 2R0/Re); whereas it enters only 
indirectly as the less sensitive abscissa for direct-force determinations 
of local friction. For example, a 10-percent error in Re easily may 
exist for most measurements presented in figure 20, and this woxild bring 
about an 8-percent error in a momentiim-method determination of Cp/Cp.; 

whereas it would bring about only 2-percent error in a direct-force 
determination of cf/cf^. 

Direct-force measurements of skin friction, on the other hand, may 
be slightly in error if a small floating element is employed because of 
the "gap effect” which at present has not been precisely evaluated. This 
effect is indicated by the available data (refs. 15, 38^ 3-nd 48) to be 
only a few percent. The present measiirements on cylinders are believed 
to be free of any such effect. Consequently, it appears that the 
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conventional momentum method is not reliable for accurate determinations 
of friction. The exact reasons for this unreliability have yet to be 
established. 

CONCLUDING REMARKS 


From a practical viewpoint the principal results of the present 
research can be summarized concisely: (1) The skin-friction coefficient 

for cylinders in axial flow is not greatly affected by variations in 
cylinder length-diameter ratio; (2) the ratio of compressible to incom- 
pressible friction is essentially independent of Reynolds number over the 
range of present experiments; and (3) measurements of average skin fric- 
tion by the conventional momentum method for the case of zero heat trans- 
fer and zero pressure gradient do not agree well with direct-force 
measurements . 

Although considerable attention has been paid (especially in the 
Appendix) to small differences that probably are of more interest to 
research workers than to practical engineers, sight should not be lost 
of the fact that agreement between the various direct-force measurements 
is quite satisfactory, particularly as regards the effect of Mach number 
on skin friction. The small uncertainties arising in the comparison of 
direct -force skin-friction measurements from several sources have stemmed 
primarily from the difficulties always involved in accurately determining 
an effective -length Reynolds number and in accurately knowing the proper 
reference values for incompressible flow. It is believed that both dif- 
ficulties mentioned have been met satisfactorily in the present experi- 
ments on cylinders: the first, by the laborious technique of making 

boundary-layer surveys at upstream stations; and the second, by making 
measurements at subsonic speeds on the same apparatus as was used for 
supersonic speeds . 

Apart from the above observations concerning detailed results of the 
investigation, there are several remarks that should be made here regard- 
ing application of the results. It might appear that the difficulty of 
determining an effective -length Reynolds number could be avoided by employ- 
ing momentum-thickness Reynolds number (Reg) as the reference parameter. 

In this regard, however, caution should be exercised because of the doubt 
cast on the accuracy of the conventional momentum method by the data 
comparison of figure 20. Uncertainties in the measured momentum decrement 
can arise because of probe interference, particularly near the wall, and 
because of the neglect of Reynolds normal stresses in the often-used 
relations between measured momentum thickness and skin friction. 

Primarily because of custom, the effect of Mach number on skin fric- 
tion has been presented for the particiilar conditions of constant -length 
Reyaolds number wherein free -stream viscosity is employed as reference 
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viscosity. Comparison vlth zero-heat-transfer flight data should he 
viewed in this light because the ratio of wall to free -stream viscosity 
is greatly different between wind tunnel and flighty and because experi- 
ence to date has indicated the wall viscosity of a turbulent boundaiy 
layer to be much more significant than the free-stream viscosity. It 
is possible that in the prediction of flight characteristics a more 
proper yet equally convenient reference viscosity for wind-tunnel measure- 
ments would be the wall viscosity multiplied by the ratio - that would 
exist in flight - of free-stream to wall viscosity. On this basis^ the 
expected effect of Mach number in reducing friction coefficients for 
flight conditions of zero heat transfer would be essentially the same up 
to Mach numbers near 2 ^ as is indicated by figure I9 of this report, but 
would be about a 7 -percent -greater reduction at a Mach number of h, and 
even greater reductions at higher Mach numbers. Because of this effect, 
and the effect that might be expected due to significant differences in 
the ratio of specific heats between flight and wind-tunnel-temperature 
conditions, the direct application of wind-tunnel results to flight con- 
ditions, even for the special case considered of zero heat transfer, 
should be viewed with reservation at Mach numbers near or above about if. 
Future research could be directed profitably toward finding methods by 
which wind-tunnel measurements of boundary layers at such Mach numbers 
can be utilized with confidence to predict flight characteristics. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
Dec. 9. 1953 
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APPENDIX 

NOTES ON DIFFICULTIES INVOLVED IN COMPARING AND ANALYZING 
VARIOUS MEASUREMENTS OF SKIN FRICTION 


There are two fundamental difficulties always involved in analyzing 
various high-speed measurements of skin friction (and, presumably^ heat 
transfer) . The first concerns the reliable determination of an effective 
Reynolds number for completely turbulent flow; such determination is 
necessary for meaningful data since the location of transition differs 
greatly between various wind tunnels as well as between wind-tunnel and 
flight conditions. The second difficulty concerns the proper reference 
values for friction in an incompressible flow; such uncertainty is becom- 
ing conspicuous because the differences in proposed formulas for friction 
at low speeds are greater than the accuracy or scatter of either Coles* 
measurements or the present measurements. Accordingly^ some discussion 
of these troublesome items is in order. 

A good illustration of both difficulties may be obtained, for example, 
by attempting to compare the effect of Reynolds number on the ratio of 
compressible to incompressible friction as determined from the present 
measurements and from Coles* measurements. The maximum range of Re 
obtainable in the present tests represents a variation of about 4:1, over 
which no systematic variation of Cp/Cpj^ with Re could be detected. The 
measurements of Coles represent a variation in Re of about 1^:1, though 
at somewhat lower Reynolds numbers. The combined sets of data for 
Moo = 2.50, shown in figure 21(a),® cover a variation in Re of about 
50;1. Over this razige Coles' data appear to indicate a small decrease 
in the ratio of compressible to Incompressible friction as Re is 
increased. As is evident from figure 21(b) , however, this small variation 
is less than the certainty with which the true reference values for Cf 
are now known, since the use of a different reference *^f such as that 

proposed by Schultz-Grunow (ref. 7)^ Schlichting (ref. 57)? ot Coles 
(ref. 15), yields results appreciably different from results based on the 
Karman-Schoenherr equation. It is difficult to determine which reference 
c^, is most appropriate to use. From the fundamental concepts involved, 

the soundest equation is the one used by Coles, which might be called ^the 
"complete" K^man equation inasmuch as it was first developed by Karman 
(ref. 13) but differs from the simplified form usually associated with 
his name. This complete equation is based solely on the two empirical 
veXocity-prof ile similarity laws referred to in the Introduction. The 

®Coles' data are plotted 1.3 percent higher than the corresponding tabu- 
lated values of reference 15 in order to convert his data, which were 
taken at Moo = 2.57, to the common value 2.50 used in this figure. 
Similarly, corrections of less than 1 percent have been applied to the 
present data to convert them to M„ = 2.50» 
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difficulty with this equation concerns the method of evaluating the 
four constants which appear in it. Coles made an elaborate evaluation 
of these consteints from velocity -profile measurements only. Before Coles* 
evaluation was available, the writers determined the same four constants 
from profile measurements of Schultz-Grunow (ref. 7) and Klebanoff -Diehl 
(ref. 58 )^ obtaining significantly different values. It would appear 
better to use shear measurements together with velocity -profile measure- 
ments in evaluating such constants . This apparently has not yet been 
done. As a result, the equation proposed by Coles does not agree as well 
with the available direct-force measurements of c^^^ as does, for example, 
the equation of Karman-Schoenherr . (A plot of the KdCrm^n-Schoenherr 
equation for c^^ on figure 1 of Coles* thesis will verify this.) On the 

other hand, the present measurements of average friction on cylinders 
extrapolated to two-dimensional flow (l/d = 0) yield values perhaps a 
little closer to Coles* Cp^ than to the Karman-Schoenherr Cfi- Also^ 

the recent Cp. measurements of Hughes (ref. 3) are closer to Coles’ 
evaluation of than to the Karman-Schoenherr equation. All such 

differences are only the order of ±5 percent from a mean^ but this is 
considerably greater than the apparent error of Coles’ measurements. It 
is believed that a definitive evaluation of c^^ and Cp^ for flat plates 
has yet to be made. 

Although the uncertainty of about ±5 percent in the proper reference 
cf^ for flat-plate data presumably could be removed by some definitive 
measurements of skin friction in low-speed flow, a second obstacle remains 
to a clear-cut comparison with Coles’ data that may not be removed by 
future measurements: The effective Reynolds number determination for 

Coles’ data is less certain than the apparent accuracy of his measurements 
would warrant. The method Coles used in determining the abscissa of the 
data of figure 21(a) was to employ the observed position of maximum 
shear stress, which is located near the end of the transition region, as 
the virtual origin of turbulence. Inasmuch as the data points not only 
are remarkably consistent, but represent three plate locations, two types 
of transition, and many tunnel pressiires, considerable credence must be 
attached to this method of determining an effective -length Reynolds 
number. However, a second method employed (and preferred) by Coles to 
determine the virtual origin yielded substantially different resxilts, 
corresponding to about 6-percent -greater skin friction at = 2.57^ 
and less at higher Mach numbers. This latter method involved relating 
the Cf (R q) curves of compressible and incompressible flow in a manner 
(details are given in Coles’ thesis) such that an effective -length 
Reynolds number could be deduced. A third method, not employed by Coles, 
would be to utilize the measured shear on the upstream floating element 
of his plate at a given tunnel pressure to calculate the virtual origin, 
and then to apply this to determine the value of Re to be assigned to 
shear measurements on downstream elements at the same pressure. In this 
latter method an assumption, such as the independence of Cf/cf^ as 
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function of Reynolds number, would be required. Data points of Coles for 
which all three of the above-mentioned methods can be applied are presented 
in table IV. The average of the three determinations of effective Reynolds 
number given in the last column has been employed for comparisons with 
Coles » data (except in fig. 21) since, in the opinion of the present 
writers, there is no decisive choice among the three methods. It is unfor- 
tunate that boundary -layer surveys were not made at upstream stations along 
Coles* plate, for such surveys would make comparison with the present data 
and other data more certain. 
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TABLE I.- LOCATION OF SURVEY STATIONS 


Survey station 

Distance from cone apex^ 
in. 

0 

3.11 

1 

3.1+8 

2 

3.96 

3 

U.82 

1+ 

6.07 

5 

8.09 

6 

9.93 

7 

1I+.90 


TABLE II.- TWENTY -FOUR POINT APPROXIMATION TO 
PITOT-PRESSURE EQUATION 


P 

Pt'-P 

|pu= 

p 

Pt' 

H| CVJ 

p 

Pt’ 

Pt'-P 

1.0000 

1.0000 

0.7000 

1.1412 

0.4000 

1.4131 

.9600 

1.0146 

.6700 

1.1606 

.3500 

1.4696 

.9100 

1.0344 

.6400 

1.1818 

.2500 

1.5791 

.8700 

1.0515 

.6100 

1.2 o 4 o 

.1900 

,1.6438 

.8300 

1.0700 

.5800 

1.2285 

.1200 

1.7174 

.7900 

1.0901 

.5500 

1.2546 

.0900 

1.7485 

.7600 

1.1060 

.5200 

1.2831 

.0600 

1.7791 

.7300 

1.1230 

.4900 

1.3143 

0 

1.8400 
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TABLE III.- VALUES OF 05-/(017.) FROM EXPERIMENTS ON OYLINDERS 

1 exp 

(INDEPENDENT OF L/d FOR RANGE INVESTIGATED) 


(a) From data of figure 17(a) 


(b) From data of figiare 17(b) 


M 

Gf/(Gp . ) 

•^1 exp 

0.51 

0.985 

.81 

.929 

1.99 

.745 

2.V9 

.671 

2.95 

.623 

3.36 

.578 

3.60 

.551 


M 

Cp/ (Op. ) 

•^1 exp 

0.51 

0.994 

.81 

.924 

1.99 

.757 

2.49 

.672 

2.95 

.630 

3.36 

.571 

3.60 

.552 


TABLE IV.- EFFEGTIVE REYNOLDS NUMBER (IN MILLIONS) 
FOR EXPERIMENTS OF DOLES 


Mob 

Cf 

Re 

Uncorrected 

Ugx/ V 0 

Gorrected 

(a) 

Average 

Re 

Re^ 

Re^ 

Rea 

2.57 

0.00181 

6,600 

4.84 

4.2 

136.1 

4.4 

(c) 

2.58 

.00166 

10,200 

8.32 

7.7 

10.3 

7.8 

8.6 

3.70 

.00162 

4,100 

3.54 

t> 2.6 

4.0 

(d) 

(c) 

3.70 

.00138 

7,560 

7.25 

6.4 

8.6 

7.6 

7.5 

4.51 

.00148 

3,470 

3.52 

131.8 

3.6 

(d) 

(c) 

4.55 

.00122 

6,590 

6.83 

5.6 

8.2 

6.0 

6.6 

4.50 

.00155 

2,900 

3.37 

(e) 

2.8 

(d) 

(c) 

4.54 

.00126 

5 , 24 o 

6.91 

(e) 

6.2 

6.1 

(c) 


determined from position of maximum shear (interpolated 
from ref. I5) 

Re2 determined from Rg and Cf^(Rg^) relationship (ref. I5) 


R03 determined from shear on upstream element, assuming 
independent of Reynolds number 

^Oorrection exceeds 25 percent of Uex/ve 

^Average of three values not possible or considered unreliable 
because of large magnitude of correction 
"^low not fully turbulent over upstream element at this pressure 
®Rei not determined by Doles 
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TABLE V.- VALUES OF Re FOR DATA OF FIGURE 20 


Experiment 

Re X 10“® 

Method of determining Re 

Wilson 

10 

extrapolate 6 to zero 

Spivack 

12.5 

extrapolate 6 to zero 

Rubes in-Maydew-Varga 

7 

ratio Cf/Cp taken as 
average of various theo- 
retical predictions 

Monaghan - J ohns on 

3 

extrapolate 6 to zero 

Br ini ch -D i ac oni s 

3 to l8 

extrapolate 6 to zero 

Ladenhurg-Ber shader 

2.7 

distance from end of tran- 
sition as indicated by 
schlieren 

Coles 

8.6 at Moo = 2.58 

7.5 at Moo “ 3.70 

6.6 at Moo = ^.55 

average of three methods, 
see table TV 

Chapman -Ke s ter 

10 at Moo = 3 . 8 , 

3 to 10 at Mn, = 2.0 

extrapolate 0 to zero 


SKIN FRICTION 

INCOMPRESSIBLE SKIN FRICTION 



VAN DRIEST(E) 
LI-NAGAMATSU(a»|) 
LCLEMMOW (E) 

MODIFIED FRANKL-VOISHEL 
RUBESIN 

FERRARI(E) 

LI-NAGAMATSU (aU) 
LCLEMMOW (m) ^ 

SMITH -HARROP 
ECKERT 

VAN DRIEST (I) 
LI-NAGAMATSU (a=0) 
CLEMMOW(I) 

WILSON 

MONAGHAN 

EXTENDED FRANKL-VOISHEL 
COPE 

TUCKER (E) 

YOUNG -JANSSEN 
VON KARMAN 
TUCKER (I) 


2 3 
MACH NUMBER, 

Figure 1.- Results of various analyses of turbulent skin friction for compressible flow over an 

insulated flat plate. 
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BOUNDARY LAYER 




FIRST DRAG MEASUREMENT 


SHROUD 


-WIRE TRIP 


SECOND DRAG MEASUREMENT 






n 




PLANE OF 
SURVEY—-. 
(STATION 0) 



PROBE 


-EFFECTIVE STARTING POSITION 

lULENT BOUNDARY LAYER 

f ^ 



Figure 2 .- Sketch illustrating method of determining friction force and 
effective starting position of turbulent bomdary layer. 




NACA TN 3097 


37 



A -15964 



A -15 965 

Figure 3-- Cone-cylinder models mounted in wind tunnel. 
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tutCTRICALLY 

NONCONDUCTING 

SHROUD 


Figure 4.- Apparatus used for survey across boundary layers. 
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PROBE 

NO. 

DIMENSIONS IN INCHES 

USED IN 
TUNNEL 

USED FOR DATA 
OF FIGURES 

h, 

h2 

W 


^2 

^3 

1 

0.001 

0.003 

0.016 

0.10 

0.21 

0.97 

NO. 1 

8, 12a, 18 

2 

.004 

.01 1 

.020 

.18 

.36 

.97 

NO. 2 

9a, 12a, 14 


Figure 5*- Probes employed for boundary-layer surveys. 



TRIP NO. 

d 

r 

1 

0.187 IN. 

0.010 IN. 

2 

.187 

.005 

3 

.250 

.010 

4 

.250 

.005 

5 

.500 

.010 

6 

.750 

.010 


Figure 6.- Dimensions of boundary- layer trips. 
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Figure 7 »~ Effect of various trips on transition as determined by near- 
surface pitot-pressure measxirements; M„,=2.9, Re=0.5xl0® per inch; 
tunnel No. 2. 
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(a) Dimensional profiles. 
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(b) Dimensionless profiles. 


Figure 8.- Effect of trip on velocity profile at station 0 (beginning 
of cylindrical surface over which skin friction was measured) ; 
1^=2. 0, Re=0 . 64x10® per inch; tunnel No. 1. 
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(b) Effect on skin friction; l/D = 13. 

Figure 9-- Measurements with various trips alone and in combination; 

Moo=2.0; tunnel No. 2. 
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Figure 10.- Measurements of skin friction at subsonic velocity using 
different nose shapes; M,*,=0.8l, L/D = 8. 
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(b) Skin friction; M„=2.95, Re=9.5xlO®, L/D = 13- 


Figure 11.- Attainment of near-equililrium conditions in tunnel No. 2. 
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(a) Comparison of velocity profiles; Re=0. 66x10® per inch. 
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(b) Comparison of skin friction; L/D = 13. 
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Figure 12. • 


Comparison of measurements in the two wind tunnels employed* 

H»=1.98. 
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(b) Skin friction for distributions A and B. 


Figure 13.- Absence of significant effect of a moderate variation 
in pressure distribution on average skin-friction coefficient; 
L/D = 13i t\innel No. 1. 



Figirre l4.- Development of boundary- layer profile along cylindrical 
surface; IVto=2.0, Re=0. 66x10® per inch; tunnel No. 2. 


SKIN-FRICTION COEFFICIENT, 



Figure 15 •- Effect of Reynolds number, Mach number, and length-diameter ratio on skin friction; 

tunnel No. 2. 
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(b) H„=2.0, l/d = 13 


(c) M„=2.9, l/d = 13 
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REYNOLDS NUMBER, Re 

(d) M„=2.5, l/d = 23 

Figure l6.- Absence of significant effect of Reynolds number on ^ for moderate range 

of Reynolds number. 
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(a) Referred to Karman-Schoenherr Cp^. (Plotted points represent individual measurements at 

Re nearest designated value.) 

Figure I7.- Effect of Mach number on skin friction; tunnel No. 2^ 
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(b) Referred to Keurman-Schoenherr (Plotted points represent values from curves faired 

through all measurements for a fixed and L/d.) 



Figure 17 •- Continued. 
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(c) Data of figure 11 {^) referred to experimental 
Figure IT.- Concluded. 
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Figure l8.- Comparison of velocity profiles on cylinder and flat plate; M^=2.0, Re-6xl0®. 
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INCOMPRESSIBLE SKIN FRICTION 



vn 


Figiire I 9 .- Effect of Mach number on turbulent skin friction obtained by direct-force method 

for constant pressure and no heat transfer. 
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Figure 20.- Comparison of average skin-friction determinations from momentimi method with those 

obtained from direct-force method. 
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(a) Comparison of Coles' data with present experiments. 
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(b) Coles' data referred to various c-p. . 
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Figure 21.- Effect of Reynolds number on ratio of compressible to incompressible skin friction* 

M«p2.50. ’ 
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